Giant Magnon Gap in Bilayer Iridate Sr 3 Ir 2 07 : Enhanced Pseudo-dipolar 
Interactions Near the Mott Transition 
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Using resonant inelastic x-ray scattering, we observe in the bilayer iridate Sr3lr2 07, a spin-orbit 
coupling driven magnetic insulator with a small charge gap, a magnon gap of «92 meV for both 
acoustic and optical branches. This exceptionally large magnon gap exceeds the total magnon 
bandwidth of «70 meV and implies a marked departure from the Heisenberg model, in stark contrast 
to the case of the single-layer iridate Sr2lrC>4. Analyzing the origin of these observations, we find that 
the giant magnon gap results from bond-directional pseudo-dipolar interactions that are strongly 
enhanced near the metal-insulator transition boundary. This suggests that novel magnetism, such 
as that inspired by the Kitaev model built on the pseudo-dipolar interactions, may emerge in small 
charge-gap iridates. 



PACS numbers: 74.10.-fv, 75.30.Ds, 78.70.Ck 

Identifying the hierarchy of energy scales associated 
with multiple interacting degrees of freedom is the start- 
ing point for understanding the physical properties of 
transition-metal oxides (TMOs). In most TMOs, the 
largest energy scale is the Coulomb interaction U, which 
suppresses charge motion in Mott insulators and allows 
description of the low-energy physics in terms of the re- 
maining spin and orbital degrees of freedom. In 5d irid- 
ium oxides, however, U is significantly diminished due 
to the spatially extended 5d orbitals, and the correlated 
insulating state cannot be sustained without the aid of 
large spin-orbit coupling (~0.5 eV) pQ. This additional 
interaction competes with other energy scales such as the 
crystal field and the hopping amplitude. The resulting 
charge gap is much smaller than that in a typical 3c? TMO 
or even those in most semiconductors, being on the or- 
der of 0.1 eV or even smaller [2J [3J. On the other hand, 
the energy scale of the magnetic interaction has been re- 
cently found to be of the same order of magnitude as 
that of 3d TMOs [H [5]. As a consequence, an intrigu- 
ing new hierarchy may result in which the energy scales 
for magnetic degrees of freedom surpass that for charge 
degrees of freedom, ushering in a new paradigm for the 
magnetism in 5d TMOs. 

Iridates of the Ruddlesden-Popper series 
Sr„ + ilr n 03„+i display a systematic electronic evo- 
lution as a function of the number of Ir02 layers (n) ; as 
n increases, the electronic structure progresses toward 
a metallic ground state as evidenced by the softening 
of the charge gap in S^I^Oy (n=2) and the metallic 
ground state found for SrlrOs (n=oo) [3J. The charge 
gap becomes so small already at n—2 that it cannot 
be resolved even in the optical conductivity spectrum, 
indicating proximity to the Mott transition point. Thus, 



the bilayer compound S^L^Oy provides a platform 
for investigating the nature of magnetism in the small 
U region near the metal-insulator transition (MIT) 
boundary. 

In this Letter, we report the magnetic excitation spec- 
tra of S^I^Oy measured by resonant inelastic x-ray scat- 
tering (RIXS) 6 , which show a number of features char- 
acterizing the unconventional nature of the magnetism 
in Sr 3 Ir 2 07 lying close to a Mott critical point. We ob- 
serve two anomalous features: a giant magnon gap of «92 
meV, even larger than the total magnon bandwidth «70 
meV, which demonstrates a marked departure from the 
Heisenberg model; and a very small bilayer splitting (sa5 
meV) , which is surprising in view of the "cubic" shape of 
the spin-orbit entangled wavefunction in iridates [3 [5] , 
which would suggest strong inter-layer interactions. The 
observed small bilayer splitting indicates frustration of 
the inter-layer isotropic exchange interactions. The tem- 
perature scale of the magnon gap exceeds 1000 K, indi- 
cating that the melting of the G-type collinear antifer- 
romagnetic (AF) order at ss285 K [5] is not driven by 
thermal fluctuations of magnetic moments, but rather 
by thermal charge carriers. Our analysis shows that the 
large magnon gap results from enhanced pseudo-dipolar 
(PD) interactions, which has an intriguing implication for 
the Kitaev model [TU] discussed recently in the context 
of honeycomb lattice iridates A2RO3 (A— Li or Na) in 
which the PD interactions play the major role [81 ITTUI6] . 

Experiments were carried out at the 30-ID beamlinc 
at the Advanced Photon Source. A horizontal scatter- 
ing geometry was used with 7r-polarized incident photons 
tuned to Ir L3 edge. A spherical diced Si(844) analyzer 
was used. The overall energy and momentum resolution 
of the RIXS spectrometer was about 30 meV and ±0.032 
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FIG. 1. (a) Image and (b-d) stack plot of the RIXS data recorded with q a b along high symmetry lines and q c fixed at J. 
Brillouin zone of the undistorted tetragonal unit cell (black square) and the magnetic cell (blue square) is shown with the 
notation following the convention for the tetragonal unit cell, as, for instance, in La2CuC>4. (e) RIXS spectra measured at four 
different q c 's with q a b fixed at (7r,7r). The q c 's of 0, it/3, 2-k/3, and ir correspond approximately to Z=25.65, 26.5, 27.35, and 
28.25, respectively. The red lines indicate the approximate peak positions at q c =0 and n. 



A -1 , respectively. By using a high resolution monochro- 
mator and improving the quality of the analyzer, the en- 
ergy resolution has been improved by more than a factor 
of four since we recently reported RIXS measurement on 
the single layer iridate Sr 2 Ir04 [I]. 

Figure 1(a) and l(b)-(d) show the image and stack 
plots, respectively, of the RIXS spectra recorded with in- 
plane momentum transfer q a b along high symmetry lines 
and <7 C fixed at ^. Three main features in the spectra 
are: (A) elastic/quasi-elastic peaks near the zero energy, 
(B) an intense and dispersive band in the range 90 — 160 
meV, and (C) a rather weak and broad feature above 
the dipersive band suggestive of two-magnon states. By 
comparing to the recent RIXS data [4] and theoretical 
study [17] on the single layer Sr 2 Ir04, it is evident that 
feature B with similar overall energy scale and momen- 
tum dependence of the intensity (see also Fig. 2) is a sin- 
gle magnon excitation. With this assignment, however, 



two anomalies are apparent: the acoustic branch appears 
to be absent, and the magnon gap is unusually large even 
for an optical mode. Typically, two branches of magnetic 
modes, acoustic and optical, are observed in other bilayer 
systems such as bilayer manganites Q2] and cuprates [TH] • 
To identify the acoustic branch, we scanned along the 
q c direction fixing q a b—{^^) where the maximal bilayer 
splitting is expected, as shown in Fig. 1(e). At q c =ir (0), 
only the acoustic (optical) branch has finite intensity, in 
accordance with the intuitive notion that magnons em- 
anate from magnetic Bragg spots. At intermediate q ci 
the spectrum is contributed to by both branches with a 
gradual shift in spectral weight from one branch to the 
other. We see an upward shift of ss5 meV as q c is varied 
from 7r to 0. Thus, it is seen that the acoustic branch also 
has a large gap and is nearly degenerate with the opti- 
cal branch. The small splitting of the two branches im- 
plies strongly frustrated inter-layer interactions, which, 



at first sight, seems inconsistent with the observed spin- 
flop transition driven by the inter-layer interactions [9]. 

To unravel this paradox, we first determine the origin 
of the anomalously large magnon gap. Such a large gap 
signals a marked departure from the Heisenberg model 
and that the magnetism in bilayer S^I^Oy is there- 
fore qualitatively different from its single layer variant 
Si'2lr04. A recent resonant x-ray diffraction study [9] 
establishes that Sr 3 Ir 2 07 has a c-axis collinear struc- 
ture, unlike the single layer Sr 2 Ir04 with in-plane canted 
moments [7]. These different magnetic anisotropies in 
Sr 2 Ir04 and S^L^Oy were captured in a magnetic ex- 
change Hamiltonian derived from microscopic interac- 
tions, which we will use here as well, adding to it longer- 
range interaction terms, which were also needed in the 
single layer Sr 2 Ir04 to quantitatively account for the 
magnon dispersion TJ. The resulting model contains 
intra-layer and inter-layer interactions; the intra-layer in- 
teractions read 



H ab = [js t s 3 + ts?S] + D{S^S y J - SJS*) 

(id) 

+ J2 Jz^Sj + J ^i$ > (i) 

((id)) (((id))) 

where and denote first, second 

and third nearest neighbors within each plane, and J, 
J 2 and J3 represent the corresponding isotropic coupling 
constants [see Fig. 3(a)]. The anisotropic coupling T in- 
cludes PD terms driven by Hund's exchange and those 
due to staggered rotations of octahedra [8]. The lat- 
ter also induce a Dzyaloshinsky-Moriya (DM) interac- 
tion, with the corresponding coupling constant D. Anal- 
ogously, the inter-layer interactions read 



(a) 



H c . = 



+ D c(sfsf +Z - s?sf +z )\ + J2 C SiS J+z , (2) 

(id) 

where the first sum runs over all sites in one plane, and 
the second over all next-nearest-neighbor pairs across the 
planes [see Fig. 3(a)]. Inter-layer interactions J c , r c , and 
D c for nearest-neighbors along c are complemented by an 
inter-layer next-nearest-neighbor coupling J 2c . 

In this model, the magnon dispersions are given by 



W± (q) = S*/,4|(q)-X|(q)-Y ± 2 (q) 



(3) 



with 



A± (q) = 4( J + T) + ( J c + T c ) - 4 J 2 (1 - cos q x cos q y ) 

-4J 3 (l- 72q )-4J 2c (lT7q) , (4) 
X ± (q) = 4 J 7q ± J c , Y± (q) = 4£ 7q ± D c , (5) 
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FIG. 2. Momentum dependence of the (a) position and (b) 
intensity of the magnon peak extracted by fitting the RIXS 
spectra (dots with error bar), overlaid with the fit from the- 
ory model (red solid lines). Because the acoustic and optical 
branches are not resolved in the experiment, the intensities 
for the two branches are summed in the theoretical calcula- 
tions with appropriate weights (See Supplementary Material 
for details) to compare to the experiment. 



where the upper (lower) sign refers to optical (acoustic) 
branches, and 7q = ^(cosq x + cosq y ). 

Following Ref . [5] , one can express all of the coupling 
constants for the isotropic and anisotropic exhange in- 
teractions, except the ones for the longer-range interac- 
tions (J 2 ,«/3 and </ 2c ), in terms of the three microscopic 
parameters ry, 8, and a, and thus fit the magnon spec- 
trum using these parameters. Here, ri=Jjj/U is the ra- 
tio between Hund's coupling and Coulomb correlation, 9 
parametrizes the degree of tetragonal distortion, and a 
is the octahedra rotation angle. Note that the PD inter- 
actions in the strong SOC limit are scaled by rj [8], and 
describing the deviation from the cubic wavefunction 
(J eff =l/2) and a are directly relevant to hopping ampli- 
tudes and therefore the superexchange interactions. The 
expressions for the dependences of the coupling constants 
on these parameters are provided in the Supplementary 
Material. This approach greatly reduces the number of 
adjustable parameters. Furthermore, experimentally a is 
determined to be «12° [501 HI]> and 6 is constrained to 
be in the range between 0.127T and 0.267T [S]. 

Figure 2 shows the fit of the experimental dispersion 
and intensity using the above model. We find that the 
large gap can only be reproduced when both rj and 9 are 
large, about 0.24 and 0.267T, respectively. Physically, the 
large rj can be understood as arising from the screening 
of U as the system approaches the borderline of an MIT, 
which is evident from the optical data showing soften- 
ing of the charge gap It is well known that while 
U is screened in the solid, Jh is not [2"2l |2"3"] , so that 
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TABLE I. Coupling constants (in units of meV) determined 
from fits to the experimental magnon dispersion. 

the metallic screening results in enhanced rj. The need 
for the large 9 can be seen from the fact that when 9 
is larger than (9 C (~ 0.257r), both in-plane (T) and out-of- 
plane (T c ) PD terms favor the c-axis moment [5] and thus 
there is a strong preference for the c-axis moment. This 
pronounced magnetic anisotropy is amplified by the large 
7], which leads to the sizable gap. While DM terms also 
contribute to the stabilization of the c-easy axis struc- 
ture [3], their effects in the magnon dispersion are much 
smaller than those from the PD terms. With this we find 
the resulting magnetic exchange interactions shown in 
Table 1. In addition to very large anisotropic couplings, 
we find that the nearest-neighbor Heisenberg coupling J 
is also enhanced compared to its measured [J and theo- 
retically estimated [2H[5S] values in Sr 2 Ir04. 

We now return to the discussion of the small bilayer 
splitting. Among the long-range interaction terms that 
were included (J 2 , J3, and J 2c ), J2c plays a critical role 
in determining the bilayer splitting. In a pure Heisenberg 
model, a small bilayer splitting would imply a small en- 
ergy difference between the two magnetic configurations 
shown in Fig. 3, and the coupling constants in Tab. 1 
corroborate this with J c w 4J 2c ( J 2c couples to four sites 
and thus cancels a four times stronger J c ). However, 
Eqs. (4) and (5) show that the PD interactions do not 
contribute to the bilayer splitting but considerably lower 
the energy of the G-type AF order in Fig. 3(a). 

The interesting situation arises that even if the inter- 
layer isotropic exchange interactions are almost com- 
pletely frustrated, the layers are still strongly coupled by 
the inter-layer PD and DM interactions that are respon- 
sible for the magnetic anisotropy. The strong inter-layer 
PD interactions manifested by the large gap do not con- 
flict with the small bilayer splitting since PD terms, hav- 
ing the Ising form, do not propagate magnons between 
the planes. 

The observed large magnon gap has two important im- 
plications. First, it raises the question as to how the 
magnetic order melts at a temperature scale (Tn~285 K) 
much smaller than the magnon gap (A m >1000 K). The 
rapid drop in the electrical resistivity when the system is 
heated through Tn [21 j suggests that the transport prop- 
erties are correlated with the magnetic order. However, 
the observed large magnon gap can hardly be reconciled 
with the standard single band spin-density wave picture 
with isotropic spin dynamics. Given that the charge gap 
(even unresolved in the optical data) might be very small, 
it is likely that AF order is destroyed by thermally acti- 
vated charge carriers that form magnetic polarons, whose 
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FIG. 3. (a) Ground state magnetic structure of SlOns. (b) A 
magnetic state with nearly the same energy when only Heisen- 
berg couplings are considered. 



motion is known to be particulary detrimental for an 
Ising-type magnetic order with large magnon gap that 
prevents a coherent charge propagation [26 . Whether 
this thermal-carrier-driven magnetic transition is a spe- 
cial case for Sr3lr2 07 or can be generally applied to other 
5c? TMOs with small charge gap [27J 12H] remains to be 
explored both experimentally and theoretically. 

Second, the enhanced PD interactions suggest a direc- 
tion for realization of the Kitaev model discussed in re- 
lated iridates A 2 Ir03 (A=Li or Na). In an ideal geometry 
where Ir ions sit on the vertices of a honeycomb lattice 
and are connected by edge-sharing oxygen octahedra, it 
has been shown that the isotropic exchange interactions 
are strongly reduced for the cubic J e fi=l/2 wavefunc- 
tion and the PD terms render a realization of the Kitaev 
model with a spin liquid ground state [8[ [11]. Experi- 
mentally, however, both Li 2 Ir03 and Na 2 Ir03 are known 
to have long-range order at Tjy~15 K [T^HIl], which sig- 
nals strong perturbation by the Heisenberg term. This 
is possibly due to less-than-ideal realization of the Ki- 
taev model in these compounds. Our study shows that 
approaching the MIT boundary in favor of large Jh/U 
may enhance the PD term and stabilize the spin-liquid 
ground state. In this regard, high pressure experiments 
on these iridates may be interesting. 

To summarize, we have revealed the unconventional 
nature of the magnetism in a spin-orbit entangled Mott 
insulator Sr3lr 2 07 lying on the verge of MIT. The system 
shows a marked departure from the Heisenberg model 
due to the strongly enhanced PD interactions. In con- 
trast to 3d oxides with small spin-orbit coupling that 
can be described by isotropic Heisenberg interactions 
with small anisotropic corrections, Sr3ir 2 07 exemplifies 
how a novel type of magnet can arise from a 5d oxide 
with strong spin-orbit coupling and a small charge gap. 
Our findings should have profound implications for other 
iridium compounds with lattice geometries in which the 
Heisenberg term is strongly suppressed. 
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